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Abstract 

Critical processes of ideal integrable models of Hele-Shaw flows are considered. A regularization 
method based on multiscaling expansions of solutions of the KdV and Toda hierarchies character- 
ized by string equations is proposed. Examples are exhibited in which the tritronquee solution of 
the Painleve-I equation turns out to provide the leading term of the regularization 
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1 Introduction 



A Hele-Shaw cell is a narrow gap between two plates filled with two fluids: say oil surrounding 
one or several bubbles of air. Several dispersionless (quasiclassical) limits of integrable systems 
have been found [l]-[6] which provide ideal models of Hele-Shaw flows in the absence of surface 
tension. Moreover, the same integrable structures also arise in random matrix models of two- 
dimensional quantum gravity [7j. Integrable systems of dispersionless type are solved by means 
of hodograph equations, so that generic initial conditions reach points of gradient catastrophe in 
a finite time. In the ideal Hele-Shaw models this feature gives rise to the cusp formation in the 
motion of the interface, while in random matrix models it manifests itself as the critical points of 
the asymptotic expansions for large matrix dimension N. In both cases a regularization mechanism 
of the underlying integrable models is required. As it is well known in random matrix theory [7j- 
|10j . the double scaling limit method provides a regularization scheme of the large A expansions 
which leads to models of two-dimensional quantum gravity. On the other hand, recent work (4J-[6J 
suggests the use of methods of asymptotic solutions of integrable systems [H]-[l3] to cure the 
singularities in Hele-Shaw flows. 

The present work is concerned with the regularization of a family of critical Hele-Shaw processes 
[I]- [6]. We mainly consider the case in which the interface develops an isolated finger which is close 
to becoming a cusp [l]-[5]. Then at a small scale the boundary of the finger tip is described by a 
curve 

Y(X) = P(X) y/X-v{xM), 

where (X,Y) are Cartesian coordinates, P{X) is a given polynomial and v(x,t\) is a particular 
solution of the dispersionless KdV (Hopf) equation 2dt x v = 3vv x . Here x ~ Q (t — t c ), where t is 
the physical time, t c is the critical time and t\ is a deformation parameter. The starting point of 
our analysis is the fact that v satisfies the dispersionless limit of the string equations [14] which 
characterize the one-matrix models of topological and two-dimensional quantum gravity |15j-|17|. 
Then, from the quasi-triviality property [18J of the KdV hierarchy, there is a unique solution u of 
the dispersionful KdV equation 

d tl u = ~(e 2 u xxx + §uu x ^, u = ^e 2k u {k \ (1) 

fc>0 

and its associated hierarchy such that v = w®. More precisely p3], the r-function associated with 
u is the large ./V limit of the Kontsevich integral over hermitian matrices 

r ^(*) : = 447^' e-diag^!,...,^), kt k = eTr(-Q) k , 



A N (Q): = J dX ex P Tr(A 2 6 - - A 3 ), B N (Q) := J dX exp Tr(A 2 9). 

In this KdV picture, singular Hele-Shaw processes correspond to critical points of the expansion of 
u, which in turn are associated to points of gradient catastrophe for . In our work we follow the 
same regularization procedure that is used in the one-matrix model description of two-dimensional 
quantum gravity [H]-[T0] : we apply a multiscaling limit method to obtain the leading term of the 
regularization of u near critical points. Then we use it to continue the Hele-Shaw flow on critical 
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regions. We notice that according to recent work |X9|, I20j . the multiscale regularization of KdV 
solutions gives a correct asymptotic approximation near the edges of the oscillatory zone (Whitham 
zone) which emerges at points of gradient catastrophe. We illustrate our strategy by regularizing 
the (2,5) critical finger studied in [5]. 

We also consider in this paper the critical processes of break-off and merging of Hele-Shaw 
bubbles [6j. The same analysis as in the critical finger case applies, but now the ruling integrable 
structure is supplied by the solution of the Toda hierarchy which underlies the large N limit of the 
partition function of the Hermitian matrix model 

t n := fdH exp (tr(£ *k H k )) . (2) 
J fc>i 

In this case the method is illustrated with an example of regularization for the merging of two 
bubbles. 

Our analysis makes use of the method developed by Takasaki and Takebe [21] for determining 
solutions of integrable systems by means of string equations (see also |22j-|24j). In the examples 
considered in both cases, KdV and Toda structures, the leading term of the regularization turns 
out to be provided by a particular solution of the Painleve-I equation (P-I) 

Wk = 6W 2 -£, (3) 

the so called tritronquee solution discovered by P. Boutroux [25 j . This is the same function which 
appears in the study of some critical processes in plasma [26] as well as in the analysis of the 
critical behavior of solutions to the focusing nonlinear Schrodinger equation [27]. As it is known, a 
different solution of P-I emerges in the random matrix models of two-dimensional quantum gravity 



2 Hele-Shaw flows and the KdV hierarchy 

In the set-up considered in [H [5] the cell is permeable to air but not oil. When air is injected the 
bubble develops a finger whose tip is pushed away and may become a cusp. By assuming that the 
finger is symmetric with respect to the X-axis and that the cusp is formed at the origin, then near 
the origin the finger turns to be described by a curve of the form 

Y(z):=( =1= ) Vz 2 -v, X = z 2 . (4) 



Here the subscript © denotes the projection of z-series on the positive powers, and tj-, k > 1 are 
deformation parameters. The function v stands for the distance between the tip and the origin 
and it is determined by imposing the asymptotic behaviour 



l+l 

Y(z) = 



^(fc + -)t fe z 2fc - 1 + ^ + 0(z- 3 ), z^oo, (5) 



2z 

k=l 



where the coefficient x is proportional to time x ~ Q (t — t c ). The resulting equation for v is the 
hodograph equation 

H(t,v) :=Y J {2k + l)t k r k (v)+x = 0, (6) 



k>l 
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where t := (x, ti, ti, ■ ■ •) with = 0, fc > Z + l, and are the coefficients of the generating function 

z sr^ r k (v) 



r :- 



\ 7. — 11 — Z 

v ^ u k>0 



-^T' r ° = L 



The hodograph equation is the basic piece to solve the system of string equations 

(z 2 )_ = 0, (mz- 1 )_=0. (7) 
for the Lax-Orlov functions of the dispersionless KP (dKP) hierarchy 

oo 

z = p + u + O(p~ 1 ), m = Y,(2k + l)t k z 2k + x + 0(z~ 2 ), (8) 

k=l 

{z, m} := z p m x - z x m p = 1. (9) 

Here and henceforth the subscripts — and + denote the projection of p-series on the strictly negative 
and positive powers, respectively. According to Theorem 1.5.1 of [21J, if (z,m) are solutions of ([7]) 
satisfying the asymptotic forms ([8]), then they verify the dKP hierarchy. More concretely, as the 
first string equation means that z 2 = p 2 + v, then the function v verifies the dispersionless KdV 
(dKdV) hierarchy dt v = 2d x rj + \. 

To solve the second string equation (mz _1 )_ = and satisfy the asymptotic behaviour (JH), one 
sets mz~ l = Y^k=i& k + 1) tk(z 2k ~ 1 )+- Thus by taking into account that (z 2k ~ 1 ) + = (z 2k ~ 2 r) 9 p, 
it follows that ([7j) reduces to ©. In particular from (j4]) and ([5]) we may identify 

Y(z) = 
v ' 2z 

so that the dynamics of the curve Y = Y(X) with respect to t is governed by the dispersionless 
KdV hierarchy. The problem is that near critical points (i c , v c ) 



dH 

dv 



gm-i H 



t c ,v c dv 



m—l 



d m H 



t c ,v c dv 



tc,Vc 



+ 0, m > 2, (10) 



the solutions v of © are multivalued and have singular derivatives (gradient catastroph). These 
situations correspond to the critical regimes of the Hele-Shaw fingers described by 

In order to find regularizations of these Hele-Shaw flows we consider solutions of the disper- 
sionful version of the KdV equation 

d h u = 7 {t 2 v>xxx + Guux). (11) 



and the higher members of its hierarchy dt i u = 2d x Rjj r \(u). Here Rj are the GeFfand-Dikii 
polynomials determined by 

e 2 (RR xx - 1 -R 2 )-2(z 2 -v)R 2 + 2z 2 = 0, i? = E^' ( 12 ) 

n>0 

or by the third-order differential equation 

d x Rn+i= (^e 2 dl + vd x + ~u x }R n . (13) 



Our first observation is that there is solvable dispersionful version of the string equations ([7|) 
given by 

(L 2 )_ = 0, i(ML-'-^- 2 )_ = 0. (14) 
Here L and M are Lax-Orlov operators ([L,M] = e) 

L = ed x + u 1 d~ 1 + ■■■ , M = ^(2j + l)tjL 2 i + t + O(L- 2 ), (15) 

of the dispersionful KP hierarchy. Now, the ± parts of a pseudo-differential operator denote the 
truncations of <9 x -series in the positive and strictly negative power terms, respectively. According 
to Proposition 1.7.11 of [21], given a solution (L,M) of (fTil) satisfying (fT5j) and [L, M) = e, then 
they are Lax-Orlov operators of the dispersionful KP hierarchy. The first string equation in (|14|) 
constitutes the KdV reduction condition and leads to a Lax operator of the form 

L = (e 2 d 2 x + u) 1 *. 

The second string equation together with the asymptotic condition on M and [L, M] = e can be 
satisfied by setting 

oo 

ML- 1 - e -L- 2 = Y,(2j + l)t J (L 2 '- 1 ) + . 

3=1 

Thus, by taking into account the identity [L 2 , (L 2j+1 ) + ] = — 2ed x Rj+i, the problem reduces to 
finding solutions of the form 

u = ^2e 2k u^ k \ (16) 



fc>0 



verifying 



or equivalently 



^2(2k + l)t k R k + x = 0, (17) 



k>l 



-^-V z R+x = ^ (18) 



where V denotes the function V := Y2j>i z 2 ^ 1 tj an d 7 is a large positively oriented closed path. 
To solve this equation for u one may use (|12p to determine the coefficients of the e-expansion 

R = y e 2n R^, RM = ? (19) 

and substitute them in (fTH|) to get the system 



V s (z)R®(z)+6iqx = 0. (20) 

From (|12p . (|19p and (|20p an iterative scheme for obtaining the expansion (|16p of it follows. 
In particular, for I = it reduces to the hodograph equation §§§ for u". This means that the 
expansion (|16p is not valid near critical points (|10p so that a different expansion must be used to 
construct a regularization of the solutions of © on critical regions. 



/ 
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3 Multiscaling expansions and asymptotic matching 



Given a m-th order critical point (t c ,v c ) (|10p of ([6]), let us introduce a new small parameter e and 
new variables tj and x given by 

e:= e 2^fT, ^ = t c j + e™f 3 , x = x c + e m x. (21) 
Let us now look for solutions to (1181) of the form 



+ ^e n u (n) , (22) 



« = V c 

n>l 

where now R is expanded in the form R = Xm>o ^" ^ • determine the coefficients of R, we 
first observe that ed x = e 1//2 cfe , so that (|12p can be rewritten as 

e (i? i?^ - - flg) - 2 (z 2 - «) i? 2 + 2 z 2 = 0. (23) 

From this equation and taking into account that v Cj % = one deduces a recursion relation for the 
coefficients R^ n > and that they can be expressed in the form 

£(«) = 5(0) y fl(0) = (24 ) 

where the functions G nir are differential polynomials in u^ k \ < k < n — r + 1. In particular, 
from (|23p it follows that the leading coefficients G n := G„ in , Go = 1 satisfy the same recursion 
relation as that arising from (|12|) (with e = 1) for the Gel'fand-Dikii polynomials R n . In other 
words 

G n = i? n (^ 1 )). 



If we now substitute (|21j) - (|22j) in (|18|) and identify coefficients of e-powers we get the system of 
equations 



r a* „ 
J ^— V ^ R[k)+6 ^ x c = 0, k = 0,...,m-l, 

i -p- v z {tc) R {k) +l-^r- v z {t) i?( fe - m ) + 5 km x = o, k 



(25) 



> m. 



Since v c is a m-th order critical point of (|1U|) we have that 

dz Vz ^ tc \k+i +hoX c = 0, k = 0, . . . ,m - 1. 

(z 2 - v c ) 2 

Hence, in view of (|24p the first m equations of the system are identically satisfied, while the 
remaining ones determine recursively the coefficients vF 1 ' for n > 1. In particular for k = m we 
get the following differential equation for the leading contribution {tW 

^ Cj m (v c ) t c A i? m (u (1) ) + ^2 c jo( v c) tj + x = 0, (26) 
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where R m (u^) is the m-th Gelfand-Dikii polynomial in uW and 



c jr (v c ) = (2j + 1) <j) — — 



dz z 2 i 



2m tJi 



From (|26p it follows that depends on the rescaled variables t through the linear combination 

J2j>i c jo(v c ) tj + x, so that (%. u (1) = c j0 (v c ) up. 

Having in mind the applications to Hele-Shaw flows, we must match the solutions (|16|) and (|22|) 
for tk = t c k, \/k > 1 and x varying on some overlap interval such us as e — ► + we have x — x c — > 
and x — ► — oo. To this end we observe that since n*- ^ satisfies the hodograph equation (j^J), then 
near an m-th order critical point (t c ,v c ) it behaves as 



u 



(0) 



(z, t c ,i,t c ,2, . . .) ~ «c + \/ c ( x ~ x c), c := -m! (-^-^-(t c ,v c )J . (27) 



Hence, the solutions it = it(°) + 0(e) and u = -u c + ew' 1 ' + 0(e 2 ) match to first order in e provided 
u^ 1 ) is a solution of the differential equation (j26[) such that 

~ V7$, x -» -oo. (28) 

Let us illustrate our analysis by studying the (2,5) critical finger considered in [5j. If we set 
tj = for all j > 2 except t 3 = 2/7, the hodograph equation © for u = becomes 



-u 3 + -hu + x = 0, (29) 



so that we have 



,(0) 



2\ 2/ V , : 7 r A 1 



5(4ii 3 + 5x 2 ) - 5x) y - r . (30) 



y/5(4ti 3 + 5x 2 ) - 5x 
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We consider the case t\ < which leads to a cusp of the type Y 2 ~ X 3 which can not be continued 
|5|. The hodograph equation (|29p has a 2-nd order critical point at 



I, .1 \/-^*c,l, V c = <J--t c>1 , 



so that the solution behaves as 



n~f c +W-— — (x-x c ), x ^ x~ 
V 15 f c 

If we take t\ = ii c (i.e. ii = 0) equation ([26]) becomes 



*g + 3^ (1)2 = -^x, (31) 



and the matching condition requires 



7(1) 



15 v c 



x. (32) 



By introducing the rescalings 



1 /5iv\ 2 / 5 ,fn / 4 \ V5 



W:=-- ^V'fiW, £:= , 



we have that the differential equation (|3ip reduces to the Painleve I (P-I) equation ([5]), while the 
matching condition reads 

W ~ £^oo. (33) 

As it is known [28], the tritronquee solution discovered by P. Boutroux [25] is the unique solution 
of P-I having no poles in the sector |arg£| < 4ir/5 for sufficiently large |£|. Moreover, it has no 
poles on the positive real axis [29] . 

One may now proceed \27\ [29] by taking a numerical approximation to the tritronquee solution 
for large positive values of £ and use it to supply initial data for the P-I equation. To this end we 
have taken the approximation provided by Eq.(3.3) of [28], and have set t\ = i C) i = — | so that 
x c = tJ = 0.64 and u c = g. Thus, with the aid of the ODE solver of Mathematica we construct 
a numerical solution that matches with the dispersionless approximation f|30j) . The numerical 
analysis shows that for e = 10 -5 we have that — (v c + eu^)\ < 5 • 10~ 4 in the x-interval 
(0.6365,0.6395). As both functions take values larger than 0.8, the relative error on this interval 
is smaller than 0.000625. The matching between both solutions on the x-intervals [0.58, 0.64022] 
and [0.6,0.6403] can be observed in Fig.l. 



Figure 1: (solid line) and u c + evP~> (dashed line) 



In this way, near the critical point the regularized Hele-Shaw dynamics associated to our asymp- 
totic approximation is given by the curve 

Y(X) = (x 2 + ^X + ^ 2 - vC^L (34) 

The function u is given by u^(x) for x < 0.638, and by v c + eu^ for x > 0.638. It is defined 
until a certain x* which corresponds to the first pole of the tritronquee solution on the negative 
£-axis. Figures 2 and 3 exhibit a sequence between x = 0.6 and x = 0.6402302. Notice that u 
is a decreasing function of x and that a cusp is formed when one of the roots of the polynomial 
P(X) = X 2 + t^X + |u 2 — | coalesces with u. Thus, for x < 0.64 the tip of the finger moves to 
the left, and it starts forming a cusp as x is close to a certain value slightly higher that 0.64 (u 
coincides with the largest root of P). Next a bubble appears at the tip of the finger. Subsequently 
a new cusp forms in this bubble (u coincides with the smallest root of P), and a second bubble 
grows from this cusp while the first bubble declines until it annihilates. Finally, the remaining 
bubble is absorbed by the finger (when both roots of P coincide). 
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Figure 2: Cusp formation and creation of a bubble which develops a new cusp 




Figure 3: Emergence, annihilation and absorbtion of bubbles 



3.1 Hele-Shaw flows and the Toda hierarchy 

In the Hele-Shaw set-up considered in [6] air is injected in two fixed points of a simply-connected 
air bubble making the bubble break into two emergent bubbles. Before the break-off the interface 
oil-air remains free of cusp-like singularities and develops a smooth neck. The reversed evolution 
describes the merging of two bubbles. 

The analysis of [6] concludes that after the break-off the local structure of a small part of the 
interface containing the tips of the bubbles falls into universal classes characterized by two even 
integers (4n,2), n > 1, and a finite number 2n of real deformation parameters t^. By assuming 
symmetry of the curve with respect to the X-axis, the general solution for the curve and the 
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potential in the (An, 2) class are 

Y(z) :=( E ^ l(fc + ' )tfc+ ^ fc ) y/(z-a){z-b), X = z. (35) 
^ y/(z—a){z-b) y © 

where a and b are the positions of the bubbles tips. Due to the physical assumptions of the problem, 
the expansion 

2n oo v 

Y(z)=J2(k + l)t k+1 z k + J2^ z->°°, (36) 

k=l k=0 

of the function Y must satisfy two conditions Yq = t (physical time) and Y\ = which determine 
the positions a, b of the tips. As it was shown in [6], imposing these two conditions on ([36]) leads 
to a pair of hodograph equations 

oo oo 

^2kt k r k ^i(u,v) = 0, ^2kt k r k (u,v) + 2x = 0, (37) 
fc=i fc=i 

where ifc = for k > 2n and 

r k (u,v) 



E 



r : 

* / ( V All 

k>0 

(38) 



a : = u — 2 y/v, b = u + 2 y^- 

These equations arise in the dispersionless AKNS hierarcy [6] . However, it is straightforward to 
see [22] that by setting 

2n 

Y = 2m-^2(k + l)t k+1 z k , x = 0, t = t h tj = 0, Vj > 2n + 2, (39) 

k=l 

they also coincide with the hodograph equations which solve the string equations 

z = z, m = m, (40) 

for the Lax-Orlov functions 



z =p + u + O{p~ 1 ), m = V kt k z k ~ l + - + 



-2\ 



Z 

k=l 

oo 

z =^ + Wl + o(p), m = Y j j i k z k ~ l - % + 0(z- 2 ), 
y =kl 

of the dispersionless 2-Toda (d2-Toda) hierarchy. In particular the first string equation represents 
the 1-Toda reduction z = z = p + u + v p^ 1 . 

In order to regularize the critical points of the hodograph equations (|37|) one may consider 
appropriate solutions of the dispersionful version of the Toda hierarchy . The natural candidates 
are provided by the string equations 

L = L, M = M, (41) 
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for the Lax-Orlov operators [L, M] = [L, M] = e of the dispersionful 2-Toda hierarchy 
The first equation determines the 1-Toda reduction L = L = A + u + v A -1 , (A := exp(ed x )). 
The system (|4ip characterizes the partition function of the hermitian matrix model in the large- N 
limit . In this way the well-known double-scaling limit method for this matrix model can be used 
to regularize the critical points of (|37p. 

As an example let us analyze the critical process of a merging of two bubbles studied in section 
VII of [61 . Thus we set ti = 0, t n = 0, n > 3, t := t\ so that (J371) reduces to 



t + 3t 3 (u 2 + 2v) = 0, 6t 3 v u + x = 0. (42) 
There is a 2nd-order critical point: v c = u 2 , with u c satisfying 

t c + 9t 3 u 2 c = 0, 6t 3 u 3 c + x c = 0, (43) 

and consequently 

4tj! + Slt 3 x 2 c = 0. 

On the other hand, the system of string equations (|4ip of the dispersionful 2-Toda reduces to 
t + 3t 3 (u 2 + v + v(x + e)) = 0, 3t 3 (u + u(x-e))v + x = 0. (44) 
where (u, v) are characterized by expansions of the form 

u = Y j e k u( k \ v = ^2e 2k vW. (45) 

fe>0 fc>0 

Obviously, from (|44p it follows that the leading terms (u^°\v^) satisfy the hodograph equations 
(J37J). To deal with solutions of (jUJ) near critical points (t c , u^°\ v^) = (t c ,u c ,v c ) of ([371) . one 
introduces 

e:=es, x = x c + e 4 x, t = t c + e 4 t. (46) 
It can be proved that (|44p admits solutions of the form 

oo oo 

u = u c + J2 £ juij) > w = ^ + S e 2j V^ 2j) . (47) 

3=2 3=1 

Thus, by equating the coefficients of e J for j = 2, 3 and 4 in (|44|) we get 

= _i F (2) ^(3) = _J_vf \ 

u c 2u c 
2(V^+u c U^) =-±- (U^f -\vtl (48) 

2 u c 6t 3 u c 

This provides us with an expression of in terms of (V^ 2 \ V^) and implies that verifies 

^ 2, + > ,2,)2 ^ (i -" A (49) 
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Near the critical point the solution of (|42p behaves as 



u ~ u c - ^Y^~(^-*) ; u ~ + -^(* c - *) as i i~ 
so that matching requires a solution of (|49p satisfying 



~ ^J--, as i--oo. (50) 
3 y *3 

Now, if we set x = x c = 1 and introduce the change of variables 

it follows that TV must satisfy the P-I equation © and the asymptotic condition W ~ ~\[\ as 
£ — ► do, so that it must be the tritronquee solution of P-I. 

Thus near the critical point the regularized Hele-Shaw dynamics of this example is characterized 
by the curve 

Y(X) = 3t 3 (X + u)y/(X -n) 2 -Av, 

where 

u = u c y (2) , v = v c + e 2 V {2 \ 

u c 

The resulting process represents the merging of the tips of two bubbles. It turns out that the right 
bubble develops a cusp, then a new bubble appears at this cusp and it grows until it merges with 
the tip of the left bubble. 
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